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Abstract. We show that the time frequency analysis of the autocorrelation 
function is, in many ways, a more appropriate tool to resolve fractional revivals 
of a wave packet than the usual time domain analysis. This advantage is crucial 
in reconstructing the initial state of the wave packet when its coherent structure 
is short-lived and decays before it is fully revived. Our calculations are based 
on the model example of fractional revivals in a Rydberg wave packet of circular 
states. We end by providing an analytical investigation which fully agrees with 
our numerical observations on the utility of time-frequency analysis in the study 
of wave packet fractional revivals. 

(Figures in this article are in colour only in the electronic version) 
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1. Introduction 

In certain quantum systems with nonlinear energy spectra, a suitably prepared wave 
packet will, in the course of its evolution, regain its initial form periodically. This is 
known as the revival of the wave packet. At some intermediate times, the evolving 
wave packet will break up into a set of replicas of its original form. This is known as 
the fractional revivals of the wave packet [1], [2] . It has been shown that the phenomena 
of revival and fractional revival occur in the wave packet dyanmics of various atomic, 
molecular and optical systems such as Rydberg atoms [31 [H 03 [7J [HI M EH ITT] , 
optical parametric oscillators [TH [T31 [TT] , the Jaynes-Cummings model [JSl [TBI H3 > 
transient signals from multilevel quantum systems [TB] , potential wells [TH1 [5D1 [5T1 [55] 
and molecular vibrational states [23l [24l [25l [26] . Extensions to systems for which the 
energy spectrum depends on two quantum numbers, have also been made in recent 
years [27j [28l [29l [30l (3TJ [32] . These phenomena have been experimentally observed 
in both atomic [6l El [101 EQ ana - molecular systems [25] . 

A widely used method for probing the revival dynamics of wave packets is based 
on a study of the autocorrelation function [4j [24J, [33] . This method is directly related 
to the observable signal in the pump-probe type experiments for studying the wave 
packet dynamics. The autocorrelation function for an evolving wave packet ^(r, t) is 
given by the overlap integral A(t) — (\&(r, t) |\&(r, 0)) . 

The auto-correlation function is a time series whose Fourier Transform (FT) will 
reveal all the frequencies, but will be unable to provide any information on when 
a particular frequency appears. This is important in the present context as the 
frequencies are time- varying. In other words, fractional revivals of a particular order 
occurs at particular instants of time. What is really desirable is a time-frequency 
analysis of the auto-correlation function such that we not only know all the frequencies, 
but also get information on when a particular frequency occurs. This is the objective 
of the present work. 

For the purpose of time- frequency analysis, short-time Fourier transform (STFT) 
has often been used in the literature. This method divides the whole time series 
in several windows, each of certain fixed width. Then the FT is performed in 
each window for obtaining the frequency information. Unfortunately, the time- 
frequency information obtained by this method has not always been satisfactory as 
its fixed window length does compromise on the frequency resolution. The method of 
continuous wavelet transform (CWT) 31 3o1 [37)ll3"7[ f3> 39, 40! overcomes the preset 
resolution problem of STFT by using a variable length window. This transform, by 
design provides good localisation in both time and frequency. The subject area of 
wavelets, developed mostly over the last fifteen years, is at the forefront of much 
current research in pure and applied mathematics, physics, computer science and 
engineering. This transform has emerged over recent years as a powerful time- 
frequency analysis. A narrow window is used for the analysis of the high frequencies 
and gives a better time resolution. A wider window is used for the analysis of low 
frequencies and gives a better frequency resolution. The continuous wavelet transform 
(CWT) of a signal f(t) is defined as 



This transformed signal is a function of two variables, s and r that are used respectively 
to scale and translate the wavelet window whereas <fi* is the complex conjugate of the 




(1) 
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transforming function known as the mother wavelet for the CWT. In our study we 
have used the Morlet wavelet as the mother wavelet. The contribution to the signal 
energy at the specific scale s and location r is given by the two dimensional wavelet 
energy density function known as the scalogram: 

E(s,t) = \T(s,t)\ 2 . (2) 

The frequencies are inversely proportional to the scale parameter, and thus s and 
r together help provide information in the time-frequency plane. In this paper, we 
present a case study to show that a wavelet based time-frequency analysis is superior 
in many ways to the standard time-domain analysis of the auto-correlation function. 



2. Fractional Revivals of a Rydberg wave packet 

We consider a Rydberg wave packet which is a superposition of circular Hydrogenic 
states having I = m = n — 1 [5] . The time dependent wave function for a localized 
wave packet formed as a one dimensional superposition of eigenstates may be written 

as 

^(r,i) = ^^„(r>- lB " f (3) 

n 

As a pre-requisite for obtaining fractional revivals, we assume that the weighting 
probabilities |c„| 2 are strongly peaked around a mean value n with a spread An = 
n max — n m i„ <C h. This allows us to expand the energy eigenvalues E n = — (2n 2 )^ 1 
(in atomic unit) in a Taylor series in n as follows. 

E n = E n + E' n {n -n) + ^(n - nf + ^E n (n - nf + .... (4) 

Neglecting the overall time dependent phase and considering up to the second 
order term, we may write E n as 

„ f (n — n) (n — n) 2 1 , , 

En = 27T ^—L + K - > , (5) 



where each term in the expansion defines important characteristic time scale that 
depend on n ; 

_ 2tt _ 2?r 

Tcl = Wy Trcv = W ( ) 

Since the energy spectrum is known in this case, one obtains 

T c i = 27m 3 , T rcv = 2T cl n/3. (7) 
The absolute square of autocorrelation function is 

f(t) = \A(t)\ 2 = y\c n \ 2 \c m \ 2 e- iE ^\ (8) 



where E nm = E n - E m . 

Fig. QJa) shows a plot of |A(£)| 2 as a function of time. This plot was generated 
by choosing |c„| 2 as a Gaussian distribution with n — 320, An = 40 and a FWHM 
given by o~ = 2.5. 

Peaks appearing in Fig.Q](a), are the signature of revivals and fractional revivals. 
Fast Fourier transform (FFT) of this time series data gives the individual spectral 
components in the frequency plane as shown in Fig. [1] (b). Although we can get 
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Figure 1: (a) Autocorrelation function of a Rydberg wave packet with n = 320, An = 40 
and a = 2.5. Time t is in a.u. (in the unit of 10 10 ). (b) Spectral components present in 
the autocorrelation function of the Rydberg wave packet. Frequency is in a.u (in the unit of 
lO" 8 ). 

complete frequency information in this way, we do not have any idea on which 
frequency appears at what time. On the other hand, the autocorrelation time series 
can be recovered by using the inverse FFT, but the frequency information goes away 
completely. 

Our goal is to acquire some frequency information in some particular times 
of interest. In the next section, we make use of continuous wavelet transform to 
investigate how one is able to resolve time and frequency in a better way. 



3. Time- Frequency Analysis 



Wavelet based time-frequency representation or scalogram of the time series data 
f(t) = |^4(£)| 2 is shown in Fig. O The scalogram was computed by using the Time- 
Frequency Tool Box (TFTB) for MatLab [41] . Here, we have used the Morlet wavelet, 
described as a complex exponential modulated by a Gaussian envelope. It is a function 
of time and given by (p(t) = Tr~ 1 l i e luJat e~ t / 2 , with the central frequency coo- Let At 
and Auj be the RMS duration and bandwidth respectively of the mother wavelet (p{t), 
where At is given by 



At = 



fZo(t-to) 2 Wt)\ 2 dt 



(9) 



The term inside the square root is the second moment of the wavelet centered at to . 
Similarly the bandwidth of the wavelet is 



Alj 



J-!L(^-^o) 2 |yHI 2 ^ 



(10) 



This mother wavelet is then used to build a set of daughter wavelets by translating 
ip(t) in time, and by dilating or contracting <p(t), which not only adjusts the mean 
frequency but also the spread of the daughter wavelet. 

Consider the case, when the mother wavelet is scaled by s. The Fourier transform 
of <p(t/s) is \s\<f(soj). The RMS duration becomes At(s) = \s\At and the corresponding 
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RMS bandwidth is Aui(s) = Aw/s. It implies At(s)Aui(s) = AtAui, which is 
independent of the scaling parameter s. Also, the RMS bandwidth is not affected 
by the translation parameter r because translating a function does not affect the 
magnitude of its FT. For our chosen mother wavelet ip(t), Eqs. © and (JTUJ) yield 
At = Auj = 1/V2 so that Am At — 1/2, which is independent of scaling and 
translation. Thus the uncertainty of time localisation is accompanied by an increase 
in the uncertainty of frequency localisation or vice versa. 

Referring to Fig. [2j we note that several patches appear in a rectangular array on 
the time-frequency plane. Each patch is centered about a particular frequency and a 
particular time. To help us understand the occurrence of these patches, we undertake 
an analytical approach as described below. Recall that the absolute square of the 
autocorrelation function is given by Eq. (|8|). Writing y = (t — r)/s, the CWT of f(t) 
can be written as 

/oo 
f{ys + TW{y)dy. (11) 
-00 

where 4>(y) is the Morlet with shifted time r and scaled by s: 

4>{y) = Tr-^V^e^ 2 / 2 . (12) 
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Figure 2: Time-frequency representation of the autocorrelation function of Rydberg atom. 
Time t (in the unit of 10 10 ) and frequency (in the unit of 10 -8 ) are in a.u. 

Substituting in Eq. (jlip and performing the integration over y, we get 

T(t,s) = v^^T \cr\ 2 \c m \ 2 ^ lj ^- lE - Tln e- { ^ sE - lhf 12 ■ (13) 

Maximum values of T(r, s), corresponding to a particular scale parameter s should 
occur whenever the factor e -(^o+sE nm /h) /2 a pp roac hes unity. This gives rise to a 
constraint: 

ujq = -sE nm /h. (14) 

Since the distribution functions |c„| 2 and |c m | 2 are peaked about n, the central 
frequency of each frequency band in Fig. [ljb) can be obtained by setting n = n 
and m = h + p, where p is an integer. Since loq is positive, we must insist that p is a 
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positive integer. Substituting in Eq. (I14[) and using the quadratic approximation ([5]) 
for the energy eigenvalues, we immediately get 

= (l-p^-) (15) 



Since the scale parameter s is related to frequency / by the relation s = Wo/(27r/) 
and T c \ <C T Tev , we finally obtain the simple formulae 

fp = 7fr, s = 7r-r- (16) 

J-cl ^Jp 

The above formula correctly predicts the central frequencies around which the spectral 
components are clustered in the frequency plane as was shown earlier in Fig.QJb). Note 
that each horizontal band in Fig. [2] has a spread about its central frequency. We will 
now show that the terms corresponding to these frequencies add up coherently in the 
expression for T(r, s) at a particular time r given by 

r = |T rev (17) 

where k is an integer. 

Coherent addition of terms would require that the phase factor exp(— iE nm T /Ti) 
be the same for these terms. That is, for arbitrary values of n and n' , one should be 
able to satisfy the condition 

TE n ,n+p = rE n ',n'+p + 27rA (18) 

where N is an integer number. Using the approximation (JSJ), it is now easy to show 
that r is indeed given by the expression (|17[) . An equivalent and simpler way of 
deriving this result is to insist that the phase factor exp(— iTE n _ n+p ) is independent 
of n. 

The expression (flTj) gives us the time instants at which |A(t)| 2 is peaked. 
Specifically, it tells us which frequency beats occur at what times. The lowest 
frequency beats f\ occur for p = 1 at times t/T ICV — 1/2, 1, 3/2, 2,.... Similarly, the 
frequency beats fi corresponding to p — 2, occur at t/T ICV = 1/4, 1/2, 3/4, 1,.... Thus 
the patches appearing in the bottom row give the transition frequency between any 
two consecutive levels (p = 1), the ones on the next row are the transition frequencies 
corresponding to p — 2 and so on. In this way, frequency bands, depicted in Fig. [2J 
trace out the signature of fractional revivals. 

Thus interestingly, in time-frequency plane, one can find both time and frequency 
information from the two time scales T c i and T rev by using the expressions (|17p and 
(|16[) . In fact, these two expressions provide the location of each patch in the time- 
frequency plane. As an example, let us consider the fourth maximum or patch in the 
third harmonic appearing in Fig. [21 Here, the third harmonic corresponds to p = 3, so 
the frequency information corresponding to this patch can be obtained from Eq. (116[) , 
which is /3 = 1.457 x 10~ 8 a.u., as shown in Fig. [2] The corresponding value for time 
can be obtained from Eq. (fT?)) . In this case, k = 4, and p = 3, so the patch will appear 
<i u i ' ' .I. » 

4. Advantages of the Time-Frequency Representation 



It is seen that the wavelet based time-frequency representation provides localisation 
in both time and frequency. Higher order fractional revivals and their localisation in 
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Figure 3: (a) Autocorrelation function of Rydberg atomic system with a — 2.5, (b), (c), (d) 
and (e) show \T(t,s)\ 2 for s = 1.96, 0.97, 0.65 and 0.48 (in unit of 10 s ) respectively. Times 
are scaled by revival time and \T(t,s)\ 2 is scaled by 10 6 . 

time are clearly manifest in the time-frequency plane. We showed how T rcv and T c \ 
are themselves sufficient to explain the time-frequency plane completely. 

Note also that the square of the autocorrelation function, plotted as a time series, 
does not resolve fractional revivals unambiguously. More precisely, the order of the 
fractional revival cannot always be determined. In contrast, the time frequency 
representation introduces a parameter p through Eqs. (fll>|) and (jTTJ) that clearly 
separates out the fractional revivals in a rectangular array on the time-frequency 
plane. In Fig. [3J we show how the corresponding s values filter out fractional revivals 
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Figure 4: (Color online) (a) A short-lived time series of autocorrelation function and (b) its 
wavelet based time-frequency representation. 



from the complicated plot of the autocorrelation function. For p = 1, one can find s 
using Eq. (JT5J) . This specific s value can filter out the signature of fractional revivals 
as shown in Fig. [3][b). This is also true for the higher order harmonics as shown in 
Fig. Etc), (d) and (e) for p = 2, p = 3 and p = 4 respectively. These specific values of 
s, filter out the signature of the corresponding higher order fractional revivals. 

The time-frequency representation can be useful for a wave packet that decays 
before its revival time. Although the short-time evolution of |A(i)| 2 can still be used 
to estimate T c \, no information can be gained about T rev if the energy spectrum of the 
system is not known. However, as long as the wave packet survives long enough for 
some patches to occur in the scalogram ( see Fig.0}, one can use Eq. (fT7|) to estimate 
T rcv . What makes this possible is the better resolution available in the time- frequency 
plane for the detection of fractional revivals. 



5. Conclusion 



In conclusion, we have made use of the continuous wavelet transform to demonstrate 
the time-frequency representation of autocorrelation function for the wave packet 
dynamics of a Rydberg wave packet. An analytical approach is provided to interpret 
the time-frequency plane and explain our numerical observations. We have shown 
that the time-frequency representation not only provides a complementary method of 
analyzing fractional revivals, it is a better tool in resolving fractional revivals. Finally, 
it is shown that the time-frequency representation may be able to extract information 
about the revival dynamics of a short-lived system even if the system decays before 
reaching its revival time. 
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